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Abstract 

By choosing an electron gas resting instead of drifting in the laboratory 
coordinate system as the initial state, the first order perturbation calculation 
of the previous paper (Phys. Stat. Sol. (b) 198, 785(1996)) is revised 
and extended to include the high order field corrections in the expression 
for the frictional forces and the energy transfer rates. The final expressions 
are formally the same as those in first order in the electric field, but the 
distribution functions of electrons appearing in them are defined by different 
expressions. The problems relative to the distribution function are discussed 
in detail and a new closed expression for the distribution function is obtained. 
The nonlinear impurity-limited resistance of a strong degenerate electron gas 
is computed numerically. The result calculated by using the new expression 
for the distribution function is quite different from that using the displaced 
Fermi function when the electric field is sufficiently high. 

Subject classification: 72.10 and 72.20 
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I. INTRODUCTION 



In the previous paper [T|], we reformulate the balance equation theory of Lei and Ting 
in ab initio manner. In this study, we neither distinguish the degree of freedom of the 
center of mass from the relative degree of freedom of the electrons nor adopt the effective 
Hamiltonian. All our calculations are carried out in the laboratory coordinate system and 
based on the original Hamiltonian {i.e. (12) to (21) of [|l|). We show that the S'-matrix is, 
in fact, also a functional of He, the electrostatic potential of the electrons. This feature is 
missing in the Lei-Ting formulation. It implies that the expressions of the frictional forces 
and the energy transfer rate all depend explicitly on the electric field. In the second paper 
0, this electric field dependence is studied perturbatively. To first order in the electric field, 
the results are given by (45), (65), (68) and (69) of 

By choosing a new initial state: both the electron gas and the phonon gas resting in 
the laboratory coordinate system, we revise and extend in the present paper the first order 
perturbation calculation given in paper [Q to high order in the electric field. If without a 
specific statement, we use in this paper the same notations as those in |]l[ and 0. 

II. FRICTIONAL FORCES AND ENERGY TRANSFER 
RATES: REVISION AND EXTENSION 

We recall that in the papers |jl| and we follow the current hot-electron transport 
theory to assume that at time t = — oo the system is composed of two independent 
equilibrium subsystems: a drifting electron gas with temperature and a phonon gas with 
temperature T. The corresponding initial density matrix is 

PO = PePph (1) 

with pe = -i-e~(-^'=~^ P~^^)/'^'= and Pph = -^e~^p'^^'^ . It is pointed out in the Appendix that 
in order to keep the consistency of the theory, one should abandon such initial state from 
beginning and adopt instead the initial state: both the electron gas and the phonon gas 



resting in the laboratory coordinate system. And so the initial density matrix ([^) should be 
replaced by 



= plPph (2) 



1 



with Pg = — e '^^^ i-iN)/Te _ Keeping this modification in mind, we revise the formulae of 0, 

-'e 



and calculate further the high order Feynman diagrams. 

In steady state, the balance equations are given by (1) to (5) of 0. As said above, the 
density matrix po in these equations should be replaced by the new one fP . Taking (67) of 
1^ into account, we rewrite these equations as follows 

eA^E + fe, + fep = 0, (3) 
eATE ■ V, + Wei + Wep = 0, (4) 

where 

fei = r7lI^(q)P/(q)(Pk+q,k+q-Pk,k) ('^^ (t)4+q, ^ (t)Ck a W ) , (5) 

^ kqcr 

fep = E ^qA(Pk+q,k+q " Pk,k) W V'q A . (t)Ck. (t) 5(t) ) , (6) 

V »/ kqA(T 

= -^E^(q)PKq)(^(k + q) -i5(k)) (5+(t)cU,,(t)Ck.(t)5(t)), (7) 
kqcr 

Wep = E AV(i?(k + q) - Eik)) {S+it)ip^xit)ci_,^Jt)cUt)Sit)) . (8) 

V K kqAa 

Here (■ • ■) denotes the statistical average over the density matrix p°, i.e. 

(■■■) = tr(---p°). (9) 

In (^), only denotes the drift velocity in the final steady state. 

As shown in [^, the frictional forces and (^, and the energy transfer rates (|^ and 
can be expressed in terms of the closed time path Green function G{k.t,k.'t'). G_ is defined 
by (17) to (20) of |^ with (■ ■ ■) redefined as (^. For example, the frictional force due to 
impurity scattering, fgj, is given by (22) of [P[, i.e. 



= -77 E (Pk'k' - Pkk) y_ — tr(LG:(kk', uj)) . (10) 
The Green function G' is calculated by the Keldysh diagrammatic technique. The Feynman 

k 

diagrams of order E and E are shown in Fig. 1 of |]^. The solid line — > — represents the 
electron propagator 

G°(kc.) = (11) 
V G^(ka;) y 

where 

G;;^(ko.) = ——^ , (13) 

^-E[y) - 2^ 

Go''(k^) = (Go^(k^) - G^(k^)) (1 - Vl) . (14) 
They are just (24) to (27) of 0, except that /k in (27) of is replaced by the Fermi function 

•^k = g(E(k)-/.)/Te + 1 • (1^) 

It is easy to see that this revision is the direct consequence of replacing the initial density 
matrix po by the new one p°. Repeating the calculation expounded in the revised version 
of (28), (40) and (41) of 0] are easily obtained. For convenience of latter discussion, we 
rewrite them in the form 



f(o) _ 27r 

ei 

2n 



J2nM^' - k)r(pk'k' - Pki.)L^'\k,k')5{E{k) - E{k')) , (16) 
^ kk' 



fi'^ = y E^Xk' - k)P(pk'k' - P^^)L^'\k, k')6{E{k) - E{k')) + 



kk' 



+ ^E^^l^(k' - k)r(pk'k' - Pkk)^(°nk,k')C(k,k')5'(i?(k) - E{k')) , (17) 



^kk' 



in which 



L(")(k,kO = (-er(k)E ■ Vk)"/k° - (-er(k')E • VkO"/k° , (18) 
C(k, k') = -eE • (k|r|k) + eE ■ (k'|r|k'). (19) 



The Feynman diagrams of order E'^ are given in Fig. 1. Their contributions to the 
frictional force fei can be written down by the rules expounded in 0. We obtain 

,(2) 1 ^-^ , , , , , , ,n f duj 



I](Pk'k' -Pkk)^i|^^(k'-k)|2y" ^ 

kk' cr 



X < hm 

ki -.k 



+ hm 

(ki,k2)k 

+ hm 

{k'i,k2)k' 



-eE ■ (k|r|ki)) (-eE ■ {k[\r\k')) tT[LG^{ku)G\kiu)G^{k[u)G%k'u 

-eE ■ (k|r|k2)) (-eE ■ (kajrlki)) tr(LG°(kcj)G°(k2a;)G°(kia;)G°(k'u; 
;-eE ■ (k'Jrlks)) (-eE ■ (k2|r|k')) tr(LG°(kcj)G°(k;cj)G°(k2Cj)G°(k', 



+ 



+ 



(20) 



The symbol (ki, k2)k3 means that ki and k2 not only approach each other but also approach 
to k-^. It is reminded that the highly singular matrix element (ki|r|k2) should be handled 
carefully by the trick expounded in and After a lengthy but straightforward 

calculation, we finally obtain 



fi? = y E(Pk'k' - Pkk)n.|t;(k' - k)\'L^'\k,k')5iEik) - Eik')) + 



kk' 



97r 

+ TT E(Pk'k' - Pkk)^.|^(k' - k)pL(i)(k, k')C(k, k')6'iEik) - Eik')) . 
^ kk' 

The same method can be used to calculate any perturbation term. The result is 



(21) 



2tt 



fe^ = v5:(Pk'k'-Pkk)ri.|t^(k'-k)P 

^ kk' 

xJ2-L^"~'\^,^')[C{k,k')y6^^\E{k) - E{k')) . 



(22) 



S^^^ denotes the j-th derivative of the delta function. Then the frictional force due to impurity 
scattering is J2 fei • It can be expressed in the following compact and transparent form: 

n=0 

f- = '^T.'^M^' - k)r(pk'k' - Pkk)(/k - /k')5(^(k') - E{k)) . (23) 
^ kk' 



Here, 



E{k) = E(k) - eE - (k|r|k) , 

OO 

/k = /k + E(-"^v(k)-Vk)"/k°- 



(24) 
(25) 



n=l 



where v(k) = ^eEr(k). 

In the same way, the expression for the frictional force fep and the energy transfer rates 
Wei and Wep derived in can be revised and extended to take the contribution of the high 
order Feynman diagram into account. The final results are 

= TT H l^q^l (Pk'k' - Pkk) 

kk'qA 

(k'=k+q) 

X [/k(l - /kO^qA - /k'(l - /k)(l + n^xMEik') - E{k) - n^^) , (26) 

2tt _ _ _ _ 

= TT E^^l^(k - ^')\\Eik') - ^(k))(/k - h')SiEik') - Eik)) , (27) 
^ kk' 

^ep = -^ Y.\MME{k')-E{k)) 

kk'qA 
(k'=k+q) 

X [/k(l - /kO^qA - /k'(l - /k)(l + n^xMEik') - E{k) - n^^) . (28) 

Comparing the above expressions for fgj, fep, Wei and Wep with the corresponding ones of 
1^, we find that (P^D, (^) to (^) are formally the same as those in first order in the electric 
field but the distribution functions appearing in them are defined by different expressions. 
We further note that for the special parabolic energy band, the distinction between E{k) 
and E{k) can be neglected for E{k!) —E{k) = E{k') —E{k). As a result, the energy transfer 
rate due to impurity scattering, Wei, is equal to zero and (|23|) , (|2^) to (|28| ) are identical 
to those of Lei and Ting if the concrete expression for the distribution function is not 
concerned. 



III. DISTRIBUTION FUNCTION: CONSISTENCY 
AND PARAMETERIZATION 

In the final steady state, the distribution function of the electron gas in a uniform electric 
field is given by (|25|) . 

We first notice the fact that this expression is an unsatisfactory one and requires 
improving. As an example, we consider the following problem. In weak electric field, we 
retain only the first two term of (pSj), i.e. 



/k = /k° - er(k)E . Vk/k° . (29) 

For simplicity, we consider only impurity scattering and the case where electrons are governed 
by parabolic energy dispersion. The balance equations of present paper (BEPP), i.e. 
(I), (H) and dl) to (H), are reduced to T = [| and 



eNE = -U . (30) 
Using (|29|), one calculates i^i from (|23|) and obtains 



where (A(k)) = ^-^S (^~a^k)^ -^(k)^(k) and rjr(k) is the transport relaxation time [0] 

^ ""'"^ '''■^(g)lVcig. (32) 



rtr(k) 47r/c^ 7o 

Obviously, both sides of (|30| ) are not equal. It means that expression (^) is inconsistent 
with BEPP. The origin of this inconsistency is due to the incorrectness of ignoring the 
renormalization of the electric field vertex in deriving these equations. Now, let the bare 
electric field vertex (o) in Fig. lb and c of be replaced by the renormalized one (•) shown 
in Fig. 2. We can show that to first order in electric field, BEPP is unchanged, but (pOf) is 
modified by replacing r(k) by a function Ak, i.e. 

/k = fl - eAkE ■ Vk/k° . (33) 

For the electron-impurity scattering, Ak is just rfj,(k) PJ^. Then, the consistency between 
the expression of /k and BEPP is recovered. Inspired by this success, we replace the bare 
electric field vertex in each of the Feynman diagram for G' by the renormalized one. We 
have the same balance equations but a different expression for /k which is given by 

oo 

/k = /k + E(-eAkE ■ Vk)"/k° • (34) 

n=l 

Ak is defined by an integral equation We do not follow this approach here due to the 
complexity of the mathematics involved. 



We instead follow the approach close to the current hot-electron transport theory 
to replace v(k) by a k- independent one, v, and therefore write (^51) in the form 



/k = /k°+E(-^v-Vkr/k°- (35) 

n=l 

The point of this alternate approach is to consider v as an unknown parameter and determine 
it in the following fashion: combining (p5|) with BEPP and the equations 



2^^ = ^^, (36) 

k 

2EPkk/k = iV^Vrf , (37) 

k 

and solving them self-consistently to get v, together with Tg, /i and for given A^, T and 
E. The consistency between the expression for /k and BEPP is obviously automatically 
satisfied. It implies that if v is determined in this way, (|35|) is identical with ( p4D within 
the approximation to neglect the k-dependent of A^. So, this alternative approach not only 
keep the merit of simplicity in mathematics but also take automatically the renormalization 
of the electric field vertex. 

We rewrite the series solution ( ^5|) compactly in the form 

= T^^-^^k • (38) 
1 + mv ■ Vk 

Multiplying both side by the operator 1 + mv • Vk, we get 

(1 + mv Vk)/k = /k° • (39) 

Then we obtain the following analytical expression for /k through solving the differential 
equation (p9|): 

/k = / d^e-^MO ■ (40) 
Jo 

It shows that /k is equal to a weight average of /k(0 which is defined as 

^'^^^ = exp[(Ek(0-/^)/Te] + l ' ^^^^ 



with 



Ek(e) = ^(k - mvO . (42) 

The parameter v connects with v^, the drift velocity of the electron gas in the final steady 
state, by (|37|). One can prove from ( |37D that v is just only for the special parabolic 
energy band. 



IV. CONCLUDING REMARKS 

It is remarkable that (|39|) is nothing but a simplified form of the Boltzmann equation. 
In fact, in the relaxation time approximation, the Boltzmann equation of an electron gas in 
the presence of a uniform electric field is 

^ + eE-Vk/k = -^2^. (43) 

at Ttr 

Ttr is the approximate k-independent relaxation time. In steady state, = 0. We rewrite 
in the form 



/k + CTtrE ■ Vk/k = . (44) 



It is just (^) with V = e'ETtr/m. Thus ( ^OD is a homogeneous steady solution of the 



Boltzmann equation in constant relaxation time approximation. We recall that in the 
approach of balance equation, the quantity v in (P0| ) is considered as an unknown parameter 
and determined by solving the complete set of the equations: ( ^OD with BEPP, ( |36D and (pTj). 



This trick of parameterization not only provides a method to choose a proper value of the 
relaxation time for given temperature T and electric field E, but also ensure ( ^0|) consistent 
with BEPP. 

We expand /k(0 i^i the series 

1 d^MO 



(e - 1)" • (45) 

?=1 



Substituting it into ( ^0]) and carrying out the ^-integration, we have 



/k = /k(l) + 7T 



1 d'hiO 



2 



+ 



1 5Vk(0 



1 3 



3 



+ ■■•, (46) 



in which /k(l) is just the displaced Fermi function used by the previous authors i.e. 

•^k = g(£;(k-mv)-/z)/T, _^ I ■ ^^^^ 



In view of (^) and (^21), we rewrite ( ^6]) in the form 
1 



/k fk — 2 



[mv ■ Vk^k(l)) TTTT^ + (mv ■ Vk)%(l 



OEM ^ 'J OEM 

V is the magnitude of v. It is clear that to the accuracy of first order in the electric field, /k 
coincides with Their difference is given by the right hand side of which belongs to 
the terms of high order in the electric field. It hints that the displaced Fermi function and 
hence the balance equation theory of Lei and Ting [0 is applicable only when the electric 
field is so low that the contribution from the terms on the right hand side of (^81) are small 
and can be neglected. As an example, we calculate the nonlinear impurity-limited resistance 
of an electron gas using either the displaced Fermi function p7|) or the new one (|40|). This 
problem has been studied by Lei and Ting based on the momentum balance equation 
and the distribution function (^Tf). We repeat their calculation in the strong degenerate 
hmit and the results are plotted as Ri/Rio versus v/vp curves in Fig. 3. The curve a ( 6 ) is 
calculated by using ( pOD ( ( ^Tj) ). Rio is the resistance Ri at v = 0. vp is the Fermi velocity. 
For N/V = 1 X lO^^cm"'^ and m = O.lme, vp = 7.7 x lO^m/s. The impurity scattering is 
assumed due to charged impurities with bare potential v{q) oc The Coulomb interaction 
between carriers is weakened due to the large permittivity e of the lattices and takes the 

2 

form m. e = 10. The two curves shown in Fig. 3 coincide approximately only when 

V < 0.2vp. But they are quite different when v > 0.2vp. Choosing r ~ Ips, the electric field 
E correspond to f ~ 0.2vp is estimated to be 90V/cm. 

It has been shown by Marchetti and Cai that the momentum and energy balance 
equations of Lei and Ting can be derived from the conventional Boltzmann equation (if we 
do not concern of the concrete expressions for the distribution function) . The same argument 



is applicable to BEPP if the shift of the band energy E{k.) in the electric field is taken 
into account in the collision term of the Boltzmann equation. In view of the achievement 
of the Boltzmann equation in the condensed matter physics, we believe that the approach 



based on BEPP and ( ^OD would give better predication for the hot-electron transport in the 
higher electric field. The detail of the numerical results will be reported in a separate paper. 

It is pointed out in |0 that our theory is applicable only when the electric field should 
not be too higher]^ so that the Wannier levels and the tunneling between the bands 



can be neglected. In such lower electric field, the intracoUision field effect |]T0[ can also be 
neglected. 
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^)For example, the Wannier level is observed only when the electric field is not lower than 
10^ V/cm for r ~ Ips @. 



APPENDIX: 



In this appendix, we show that there is an inconsistency in the theory presented in P|. 
The statistical average of the total momentum of electron gas P = I] PkkCkaCkcr is given 

k(7 

by 

m = T.pMctAtM)) ■ (Al) 

kcr 

Here (■ • •) = tr(- ■ ■ po). By the method of 0, we express (|A1|) in terms of the closed time 
path Green function as follows: 



(t) = -I ^ Jm ^ Pkktr (LG(kt, kt')) • (A2) 



k. 

Here, G_ is defined by (17) to (20) of 0. In steady state, G(kt,kt') depends only on k and 
t — t' and so P(t) is independent on time. Performing Fourier transformation of G with 
respect to t — t', we obtain 

P = E Pkk / |^e^-°"tr(LG(k, .;)) . (A3) 

kcr 

Calculating G_ to first order in the electric field 0, we have 



P = EPkk/|^e^^°'tr(LG°(k,c.)) + 

kcr 

+ EPkk / :^trfL limjG°(k,o;)V^(k, k')G°(k', a;)]) +0(i?2) , (A4) 
where G° is defined by (24) to (27) of O ne can easily prove the result 

P = iVmvrf + Pkk(-er(k)E) • Vk/k + OiE") , (A5) 

kcr 

in which /k is given by (12) of 0. The self-consistent condition requires that the drift 
velocity of the electron gas in the initial state, v^, should be equal to that in the final steady 
state, ]^P ||l[. We see from ( |A5|) with disappointment that this self-consistent condition 
does not hold. It implies that there are some improper things involved in the current balance 
equation theory. 



Examining the theory carefully, we find the inconsistency results from we choosing the 
drifting electron gas as the initial state. In order to keep the consistency of the theory, we 
should abandon such initial state and adopt instead the electron gas resting in the laboratory 
coordinate system as the initial state. In connection with this, the self-consistent condition 
to identify the drift velocity of the initial state with that in the final state is no longer 
required. Then the first term on the right hand side of ( [A^ ) is equal to zero, and ( [X5| ) is 
nothing but the equation defined the drift velocity of the electron gas in the final steady 
state, jj—P- The said inconsistency no longer appears. 
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FIGURES 

Fig.l Feynman diagrams for G', each with two bare electric field vertices. 

Fig. 2 (a) Feynman diagrams for G', each with one renormalized electric field vertex, 
(b) Renormalized electric field vertex in the ladder approximation. 

Fig. 3 Ri/Rio versus v/vp curves. 

Curve a, using (HO); curve b, using (H^. 
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